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Abstract 

The Day a Bay Reactor Neutrino Experiment is designed to determine precisely the neutrino mixing angle #13 with 
a sensitivity better than 0.01 in the parameter sin 2 2#i3 at the 90% confidence level. To achieve this goal, the collab- 
oration will build eight functionally identical antineutrino detectors. The first two detectors have been constructed, 
installed and commissioned in Experimental Hall 1, with steady data- taking beginning September 23, 2011. A com- 
parison of the data collected over the subsequent three months indicates that the detectors are functionally identical, 
and that detector-related systematic uncertainties exceed requirements. 

Keywords: neutrino oscillation, neutrino mixing, reactor, Daya Bay 
PACS: 14.60.Pq, 29.40.Mc, 28.50.Hw, 13.15.+g 



1. Introduction 

The precise determination of the neutrino mixing an- 
gle 613 by the Daya Bay Reactor Neutrino Experiment 
(Daya Bay) requires measurement of the antineutrino 
flux from the six nuclear reactors at different baselines 
using eight antineutrino detectors 1 1]. Detection of an- 
tineutrinos is via the inverse beta-decay (IBD) reaction: 

v e + p — > e + + n . (1) 

The positron rapidly annihilates with an electron 
(prompt signal) while the neutron first thermalizes be- 
fore being captured by a nucleus and releasing energy 
(delayed signal). 

The value of sin 2 26 13 can be determined by compar- 
ing the observed antineutrino rate and energy spectrum 
with predictions assuming oscillations. The number of 
detected antineutrinos N^ ti is given by 

Wdet = 4^2/ €<rP sm (E,L,0 l3 )SdE (2) 

where N p is the number of free protons in the target, L is 
the distance of the detector from the reactor, e is the ef- 
ficiency of detecting an antineutrino, cr is the total cross 



section of the IBD process, P^m is the v e — > v e survival 
probability that depends on the value of sin 2 26n, and S 
is the differential energy distribution of the antineutrino. 

With only one detector at a fixed baseline from a re- 
actor, according to Eq. [2l we must determine the ab- 
solute antineutrino flux from the reactor, the absolute 
cross section of the IBD reaction, and the efficiencies 
of the detector and event- selection requirements in or- 
der to measure sin 2 26^. It is a challenge to reduce the 
systematic uncertainties of such an absolute measure- 
ment to sub-percent level, especially for reactor-related 
uncertainties. 

Mikaelyan and Sinev pointed out that the systematic 
uncertainties can be greatly suppressed or totally elim- 
inated when two detectors positioned at two different 
baselines are utilized 0]. The detector closer to the re- 
actor core is primarily used to establish the flux and en- 
ergy spectrum of the antineutrinos. This relaxes the re- 
quirement of knowing the details of the fission process 
and operational conditions of the reactor. In this ap- 
proach, the value of sin 2 26^ can be measured by com- 
paring the antineutrino flux and energy distribution ob- 
served with the far detector to those of the near detector. 

According to Eq. [2] for a single reactor core and sin- 
gle near and far detectors, the ratio of the number of 
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antineutrino events with energy between E and E + dE 
detected at distance Lf (far detector) from the reactor 
core to that at a distance L n (near detector) is given by 

(3) 



EH3 



Nf 



p,n 



P sm (E,Lf, 613) 



P sm (E,L n ,6 13 ) 



where N p j and Af p?n refer to the number of target pro- 
tons at the far and near sites, respectively. The relative 
detector efficiency (ef/e n ) can be determined more pre- 
cisely than the absolute efficiency. Hence, the detector- 
related systematic uncertainty in this approach is greatly 
reduced. Furthermore, the use of multiple modules at 
each site enables internal consistency checks. Day a Bay 
will implement this strategy by deploying two function- 
ally identical modules at each of two sites near the reac- 
tor cores, and four detectors at a site further away. 

In this paper we will compare the performance of 
the first two antineutrino detectors (ADs) using three 
months of data. The detectors were installed and com- 
missioned side-by-side in the Daya Bay Experimental 
Hall (also known as EHl) during the summer of 2011. 
The data used in this analysis were collected between 
September 23 and December 23, 2011. To reduce the 
potential for biases during data analysis, Daya Bay has 
adopted a blind analysis. The baselines, the thermal 
power histories of the cores, and the target masses of the 
antineutrino detectors will be blinded. Before unwind- 
ing, nominal values for these quantities will be used. An 
overview of the experimental site and detectors will be 
provided first, followed by a detailed comparison of the 
detectors' performance. 

2. Experimental site 

The Daya Bay nuclear power complex is located on 
the southern coast of China, 55 km to the northeast of 
Hong Kong and 45 km to the east of Shenzhen. As 
shown in Fig. 1 , it consists of three nuclear power plants 
(NPPs), Daya Bay NPP, Ling Ao NPP, and Ling Ao-II 
NPP. The complex faces the sea on the southeast, and 
is adjacent to mountains on the northwest. Each NPP 
consists of two reactor cores. All six cores are func- 
tionally identical pressurized water reactors (Framatone 
M310 and its derivative CPR1000) of 2.9 GW thermal 
power yfl. The last core started commercial operation 
on Aug. 7, 2011. The distance between the cores for 
each pair is about 90 m. The Daya Bay cores are sep- 
arated from the Ling Ao cores by about 1100 m, while 
the Ling Ao-II cores are around 500 m from the Ling 
Ao cores. 

The Daya Bay experimental facility consists of sur- 
face facilities, three underground experimental halls, 
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Figure 1: Layout of the Daya Bay experiment. The dots are reactor 
cores, labeled as Dl and D2 for the Daya Bay NPP, LI and L2 for the 
Ling Ao NPP, and L3 and L4 for the Ling Ao II NPP. Two antineu- 
trino detectors, AD1 and AD2, are currently installed in the Daya Bay 
Experimental Hall (EHl). 



and two additional underground utility halls known as 
the Liquid Scintillator Hall (LS Hall) and the Water 
Hall. The surface facilities include a Surface Assem- 
bly Building (SAB) where the ADs are assembled, an 
office building, and a building housing the ventilation 
equipment for the underground halls. Each detector hall 
contains a water pool instrumented to detect Cherenkov 
radiation, either two or four ADs installed inside the wa- 
ter pool, and modules containing four layers of resistive 
plate chambers (RPCs) over the top of the pool. The 
LS Hall is the location used for producing and storing 
the liquid scintillator (LS) and gadolinium-doped liq- 
uid scintillator (Gd-LS), as well as for filling the ADs. 
The Water Hall is used to produce purified water for the 
water Cherenkov detectors. The underground halls are 
connected by horizontal tunnels with a 0.3% slope to 
facilitate drainage of water. The surface-access tunnel 
is 267 m long with a 10% downward slope. The total 
length of the tunnel system is 3100 m. The tunnels are 
6.2 m wide by 7.0 m high to allow for transportation of 
assembled ADs. 

The mountain contour over the tunnels and experi- 
mental halls was surveyed prior to the ground-breaking 
of civil construction. An additional survey was com- 
pleted later using GPS and Total Station technologies 
that determined the position of the detectors in each ex- 
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periment hall with respect to the reactor cores to a pre- 
cision of a few centimeters. Based on these surveys, 
the Day a Bay Experimental Hall (EH1) has an overbur- 
den of 280 meter-water-equivalent (mwe) and is about 
360 m from the center of the twin cores of the Day a 
Bay NPP. The Ling Ao Experimental Hall (EH2) is 300 
mwe deep and is about 500 m on average from the four 
cores of the Ling Ao NPP and Ling Ao-II NPP. The Far 
Hall (EH3) is 880 mwe deep, about 1910 m from the 
cores of the Daya Bay NPP and about 1540 m from the 
cores of the Ling Ao and the Ling Ao-II NPPs. These 
approximate values are tabulated in Table 1 . The pre- 
cise locations will be used in the final analysis. At full 
thermal power, each AD in EH1 is expected to observe 
about 800 IBD events per day where the neutron is cap- 
tured by a Gd nucleus. 





Overburden 


D. B. 


L. A. 


L. A. II 


EH1 


280 


360 


860 


1310 


EH2 


300 


1350 


480 


530 


EH3 


880 


1910 


1540 


1550 



Table 1 : Approximate values for the overburden above and distances 
to the three experimental halls from the Daya Bay, Ling Ao, and Ling 
Ao-II NPPs. The overburden is in meter- water-equivalent and the dis- 
tances are in meters. 



for example, from the photo-multiplier tubes (PMTs) or 
SSV. Three automated calibration units (ACU-A, ACU- 
B, and ACU-C) are mounted at the top of the SSV. Each 
ACU contains a LED as well as two sealed capsules 
with radioactive sources that can be lowered individu- 
ally into the Gd-LS along either the centerline or inner 
edge, or in the LS. 



ACU-B 



ACU-A 



ACU-C 




Figure 2: Schematic for a Daya Bay antineutrino detector. 



3. Detectors 

In each experimental hall there are three different 
kinds of detectors: the ADs, the water Cherenkov de- 
tectors, and the RPC detectors. In total, there will be 
eight 110-ton ADs, three water pools filled with 4400 t 
of purified water, and three arrays of RPCs covering a 
total of 800 m 2 . 

3.1. Antineutrino detectors 

Each AD has three nested cylindrical volumes sepa- 
rated by concentric acrylic vessels as shown in Fig. 2. 
The outermost vessel is constructed of stainless steel 
and is known as the SSV. The innermost volume holds 
20 1 of 0. 1 % by weight Gd-LS that serves as the antineu- 
trino target. The middle volume is called the gamma 
catcher and is filled with 21 1 of un-doped liquid scintil- 
lator (LS) for detecting gamma-rays that escape the tar- 
get volume. The gamma-catcher increases the contain- 
ment of gamma energy thus improving the energy reso- 
lution and reducing the uncertainties of the antineutrino 
detection efficiency. The outer volume contains 37 t of 
mineral oil (MO) to provide optical homogeneity and 
to shield the inner volumes from radiation originating, 



The acrylic vessel holding the Gd-LS has a diame- 
ter and height of 3.1 m and a wall thickness of 10 mm. 
This inner acrylic vessel (IAV) is nested within an outer 
acrylic vessel (OAV) with a diameter and height of 4 
m and a wall thickness of 18 mm. Both the IAV and 
OAV are made of UV- transparent acrylic @, @]. The 
lids are cone-shaped with a 3 degree tilt angle. There are 
two calibration pipes made of Teflon bellows connect- 
ing the IAV to the ACUs on the top of the AD. They are 
nested within larger Teflon bellows attached to the OAV. 
A third pipe located at the edge of the OAV provides ac- 
cess for ACU-C to the gamma catcher. Gd-LS and LS 
can flow along the calibration pipes to an overflow tank 
located at the top-center of the SSV. The overflow tank 
is a nested two-layer acrylic vessel. The mineral oil has 
two separate overflow tanks also situated on the top of 
the SSV. The SSV has a diameter and height of 5 m and 
a wall thickness of 12 mm. It is reinforced with ribs 
at the bottom, under the lid, and on the inner and outer 
surfaces of the barrel to provide sufficient mechanical 
strength for lifting after the AD is filled, to withstand 
the water pressure, and maintain an accuracy of 2 mm 
for the critical AD components and the location of cal- 
ibration sources. Two 4.5-m diameter, 2-cm thick re- 
flective panels are placed at the top and bottom of the 
OAV to increase the photon-statistics and improve the 
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uniformity of the energy response The reflectors 
are laminated with a film of ESR (Vikuiti™ Enhanced 
Specular Reflector Film) sealed between two 1-cm thick 
acrylic panels. The ESR film has a reflectivity greater 
than 98% across most of the relevant spectrum. 

There are 192 8-inch PMTs (Hamamatsu R5912) 
mounted on eight ladders installed along the circum- 
ference of the SSV and within the mineral oil volume. 
The PMT bulb is 20 cm from the OAV. To minimize 
non-uniformity of response, the PMTs are recessed in 
a 3 -mm thick black acrylic cylindrical shield located at 
the equator of the PMT bulb. A conical magnetic shield 
is wrapped around the dynode structure to protect the 
PMTs from stray magnetic fields. The photocathodes 
operate at ground potential. A single coaxial cable is 
used to supply positive high voltage and transmit the 
PMT signal to the front-end electronics. The decou- 
pling of the signal from the high voltage is performed 
inside the electronics room. The high-voltage system 
uses a CAEN SY1527LC mainframe. Each mainframe 
houses eight A1932AP 4 8 -channel high- voltage distri- 
bution modules. During data taking, the gain of the 
PMTs is set to 1 x 10 7 . 

ACU-A sits on the central axis of the detector. ACU- 
B is located at a radius of 135.00 cm to calibrate and 
study edge effects within the IAV. ACU-C is located at a 
radius of 177.25 cm for calibrating the gamma catcher, 
on the opposite side to ACU-B. Each ACU is equipped 
with a LED, a 68 Ge source, and a combined source of 
241 Am- 13 C and 60 Co. The Am- 13 C source generates 
neutrons at a rate of 0.5 Hz. The rates of the 60 Co and 
68 Ge sources are about 100 Hz and 15 Hz, respectively. 
Since the AD is fully submerged in water, the ACUs are 
operated remotely. The sources can be located to bet- 
ter than 0.5 cm along a vertical line down to the bottom 
of the acrylic vessels. When not in use, the LED and 
sources are retracted into the ACUs, that also serve as 
shielding for the sources. 

Six 2-inch PMTs (Hamamatsu R7724) are installed 
at the top and bottom of the AD, to monitor the attenu- 
ation length of the Gd-LS and LS via optical windows 
on the reflective panels. A mineral oil clarity device is 
installed on the AD lid to monitor the attenuation length 
of the mineral oil by detecting blue LED light reflected 
back from a retroreflector at the bottom of the AD. The 
AD is also instrumented with two CMOS cameras (nor- 
mally off) mounted at the top and bottom of a PMT lad- 
der, and temperature sensors along the PMT ladders and 
in the overflow tanks. The liquid levels of the LS, Gd- 
LS, and MO within the overflow tanks are each moni- 
tored using redundant sensors (capacitance, ultrasonic, 
and cameras). 



All Gd-LS ItZD and LS was produced in the LS 
Hall 1 8]. The solvent of the LS and Gd-LS, linear alkyl- 
benzene (LAB), was procured in one batch [9]. The 
fluor and wavelength- shifter are 3 g/L PPO and 15 mg/L 
bis-MSB, respectively. Gadolinium chloride reacts with 
3,5,5-trimethylhexanoic acid (TMHA) to form a solid 
complex Gd(TMHA) 3 that dissolves in LAB. The Gd- 
LS was produced in 4-ton batches and stored in five 40- 
ton acrylic storage tanks, each with internal circulation 
systems. Four tons of Gd-LS are taken from each of the 
five storage tanks and pre-mixed prior to the filling of 
each AD to ensure that all eight ADs have the same Gd 
content. The LS is also produced in 4-ton batches and 
then stored in a 200-ton storage pool with internal cir- 
culation. The MO is stored in a different 200-ton pool. 
The Gd-LS has 87.7% carbon content by weight, 12.1% 
hydrogen, and 0.103% Gd. The density of Gd-LS, LS, 
and MO is 0.860, 0.859, and 0.851 g/ml, respectively. 

To reduce systematic uncertainties and to reduce con- 
tamination from dust, all ADs are assembled in pairs in 
a large ISO 7 (class 10,000) clean room in the SAB. 
Each pair is filled over a short interval of time so that 
detector pairs are matched in Gd-LS quality and char- 
acteristics. Gd-LS, LS, and mineral oil filling occurs 
underground in the LS Hall. All three liquids are filled 
concurrently to maintain a uniform liquid level across 
all three volumes. MO and LS are pumped directly 
from the large storage reservoirs, and their mass is mea- 
sured by Coriolis flow meters. Gd-LS is first transferred 
equally from five separate Gd-LS storage tanks to a 20- 
ton ISO tank lined with Teflon and instrumented with 
load cells before being pumped into the AD. The AD 
target mass can be determined to O(0.1%) by weighing 
the ISO tank before and after the filling. All ADs are 
filled with the same filling system. 

All detector components passed low background test- 
ing. The stainless steel of the SSV was specially made 
in one batch with low-radioactivity iron ore. PMTs with 
low background glass were procured. Chemicals used 
for synthesizing the Gd-LS, GdCl 3 , TMHA, and PPO 
were purified to get rid of some radioactive contami- 
nants. The liquid scintillators are covered with a contin- 
uous flow of dry nitrogen gas to prevent oxidation while 
in storage and within an AD. 

3.2. Muon system 

The muon detector consists of a RPC tracking device 
and an active water shield. The water shield consists 
of two optically separated regions known as the inner 
(IWS) and outer (OWS) water shields. Each region op- 
erates as an independent water Cherenkov detector in- 
strumented with PMTs. The water shield as a whole has 
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multiple purposes. It detects muons that can produce 
spallation neutrons or other cosmogenic backgrounds in 
the ADs. The pool also moderates neutrons and attenu- 
ates gamma rays produced in the rock and other struc- 
tural materials in and around the experimental halls. 
The water pool is designed so that there is at least 2.5 
m of water surrounding each AD in every direction, as 
seen in Fig. [3l The near site water pools each contain 
1200 t of purified water produced at 18 MQ-cm while 
the far water pool contains 1950 1 of purified water. The 
PMTs are distributed between the inner and outer zones 
that are optically divided by Tyvek sheets. Each pool is 
outfitted with a light-tight cover with dry-nitrogen flow- 
ing underneath. 

Each water pool is covered with an array of RPC 
modules fioll . The 2 m x 2 m modules are deployed in 
an overlapping pattern (that minimizes dead areas) on a 
steel frame mounted on rails, so that the assembly can 
be retracted to provide access to the water pool. There 
are four layers of bare RPCs inside each module, with 
one layer of readout strips associated with each layer of 
bare RPCs. The strips have a "zigzag" design with an 
effective width of 25 cm, and are stacked in alternating 
orientations providing a spacial resolution of ~8 cm. 




Figure 3: Schematic for the Daya Bay Near Hall (EH1) including the 
ADs, water shields, and RPCs. 

3.3. Electronics 

Each detector unit (AD, IWS, OWS, and RPC) is 
read out by a single VME crate. There are two types 
of readout crates, PMT readout and RPC readout. All 
PMT readout crates are physically identical, differing 
only in the number of instrumented readout channels. 
Each near site has four PMT readout crates and one RPC 
readout crate, while the far site has six PMT readout 



crates and one RPC readout crate. The PMT readout is 
documented elsewhere 1 11] and is shown schematically 
in Fig.|4l In brief, the front-end electronics board (FEE) 
receives raw signals from up to sixteen PMTs, sums 
the charge among all input channels, identifies over- 
threshold channels, records timing information on over- 
threshold channels, and measures the charge of each 
over- threshold pulse. The FEE in turn sends the num- 
ber of channels over threshold and the integrated charge 
to the trigger system as well as the FADC board. The 
FADC samples and records the unshaped energy sum at 
1 GHz. When a trigger is issued, the FEE reads out the 
charge and timing information for each over-threshold 
channel, as well as the average ADC value over a 100 ns 
time- window immediately preceding the over-threshold 
condition (pre ADC). 

The RPC readout consists of 32-channel front-end 
cards (FECs) mounted on the detector modules. Each 
FEC reads out one RPC detector module. The RPC trig- 
ger module (RTM) and the RPC readout module (ROM) 
sit in a VME crate in the electronics room. Transceiver 
modules in custom crates mounted on the RPC detector 
frame relay trigger and timing information between the 
RTM and FECs, while also transmitting hit information 
to the RTM and ROM. The RPC readout is independent 
of the PMT readout and consists of a digital hit map 
of the over-threshold channels along with a GPS time- 
stamp for the trigger. Additional details can be found 
in fl2ll. 
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Figure 4: Block diagram outlining the PMT readout electronics. 

3.4. Triggers 

Triggers are primarily created internally within each 
VME crate, although the system is also designed to ac- 
cept external trigger signals. Each PMT readout crate 
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contains a local trigger board (LTB) that receives the 
number of over-threshold channels (Nhit) as well as the 
summed charge (E-Sum) from each FEE [13]. The Nhit 
and E-Sum data form the basis for two physical triggers, 
the multiplicity trigger and energy- sum trigger. In addi- 
tion to the physical triggers, the LTB can also generate 
two internal triggers, a periodic trigger and a software 
trigger. Additionally, the LTB can receive and process 
two external trigger requests, the calibration trigger re- 
quest and a cross trigger request. Cross triggers allow 
one detector subsystem to trigger other detector subsys- 
tems. Each type of trigger can be enabled or disabled 
via VME interface. The enabled triggers are ORed to 
generate a final trigger signal (local trigger) that is dis- 
tributed to the FEEs by front-panel cables. This design 
is flexible and provides various options to study and 
increase the overall trigger efficiency. The RPC sub- 
detector uses custom front-end electronics that employ 
a self-triggering scheme. A master trigger board (MTB) 
coordinates cross triggers, allowing a sub-detector to 
be triggered by another sub-detector. Calibration trig- 
ger requests are used to monitor detector performance 
and come into the system through the MTB. The MTB 
broadcasts this trigger to the appropriate LTB as a cross- 
trigger command. To avoid an occurrence of overflow 
in either the LTB data-buffer or the FEE data-buffer, the 
LTB can check if either is nearly full. If so, the lo- 
cal trigger and the corresponding data package can be 
blocked. The blocked trigger number is recorded and 
read out for calculating the dead time of the trigger- 
system offline. 

3.5. Data acquisition 

The data acquisition (DAQ) architecture is designed 
as a multi-level system using embedded Linux, ad- 
vanced commercial computer and distributed network 
technology, and is modeled after the BESIII and AT- 
LAS DAQ systems [Q. The readout DAQ operates 
concurrently within each VME crate from a MVME 
5500 single board computer running TimeSys (a real- 
time LINUX based OS). Event readout begins with an 
interrupt from the trigger module that initiates read out 
of data fragments from all modules residing in the crate 
using chained block transfer (CBLT) mode. Fragments 
are concatenated into an event ordered in time to cre- 
ate a data stream for that VME crate. The front-end 
VME system transmits the data to the back-end DAQ 
that merges and sorts the events from the various in- 
coming streams by trigger time-stamp. The DAQ sys- 
tem also provides an interactive interface for the elec- 
tronics, trigger and calibration systems. Run control is 
flexible and configurable, allowing global operation of 



all detector systems or operation of sub-sets of detectors 
whenever debugging or commissioning is required. 

3.6. Data stability and selection 

Stable data taking in EH1 began September 23, 201 1 . 
Table [2] summarizes the experimental livetime from 
September 23, 201 1 to December 23, 201 1 . The fraction 
of the physics-data-taking (non-calibration and non- 
diagnostic) time in total calendar time (2194.3 hours) 
was 86.8%. Furthermore, the data quality is good with 
1684.2 hours of data deemed suitable for physics analy- 
sis. We excluded 221 .0 hours of data from physics anal- 
ysis, including 203.2 hours of systematic studies, 10.2 
hours due to some coherent noise pickup in the elec- 
tronics, and 7.6 hours due to electronics, high voltage, 
or DAQ problems. The systematic studies primarily oc- 
curred between December 10 and December 18, 2011 
when the HV for certain PMTs were lowered to inves- 
tigate 'flasher' events (see Section 5.4). The DAQ dead 
time due to full data buffers was determined with dedi- 
cated scalers to be less than 0.0025% for both AD1 and 
AD2. 



Total calendar time (hr) 2194.3 

Total DAQ time (hr) 2092.5 

Physics DAQ time (hr) 1 905 .2 

Good run time (hr) 1684.2 

Table 2: Summary of experimental livetime. 



3. 7. Offline processing 

During data taking, the onsite experimental raw data 
are transferred to the computing centers located at 
the Institute for High Energy Physics (IHEP) in Bei- 
jing as well as Lawrence Berkeley National Laboratory 
(LBNL) in California in real time. The multi-flow con- 
current transmission mode is adopted to improve the 
data transmission bandwidth. As soon as the data reach 
the destinations, a keep-up data processing takes place. 
To help assess and monitor data quality, distributions 
of selected variables are automatically plotted and pub- 
lished on the web. After checking the data quality, the 
offline data processing reconstructs and tags the data for 
physics analysis. 

The Day a Bay offline software (known as NuWa) has 
been developed using Gaudi 1 15] as the underlying soft- 
ware framework in order to provide the full functional- 
ity required by simulation, reconstruction and physics 
analysis. NuWa employs Gaudi 's event data service as 
the data manager. Raw data, as well as other offline data 
objects, can be accessed from the Transient Event Store 
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(TES). According to the requirement from the prompt- 
delayed coincidence analysis, another specific Archive 
Event Store (AES) has also been implemented to pro- 
vide the function of look-back in time. All the data ob- 
jects in both TES and AES can be written into or read 
back from ROOT [16] files through various Gaudi con- 
verters. In this framework, reconstruction algorithms 
have been developed to construct the energy and the ver- 
tex of the antineutrino event from the charge pattern of 
the PMTs. The detector-related parameters and calibra- 
tion constants needed by the reconstruction are stored 
in an offline central database with a number of mirror 
sites located at different institutes. The algorithms can 
access the contents in the database via an interface soft- 
ware package called DBI. 



The Day a Bay simulation is based on GEANT4 111711 
with certain critical features validated against exter- 
nal data, where available, or other simulation packages 
such as MCNPX [18], FLUKA d, GCALOR 
and GEANT3. Features such as the gamma spectrum 
from a neutron capture on Gadolinium, the cosmic 
muon flux, and specific decay chains have been custom 
built. Tracking of optical photons is used during high 
precision simulations, and a comprehensive electronic 
readout simulation has been implemented. The Monte 
Carlo (MC) simulation is tuned to match observed de- 
tector response. To simulate time correlations of unre- 
lated events, different categories of simulated events are 
mixed in the output file. 

4. Low level performance 

4.1. ADC calibration 

Each readout channel is calibrated with data collected 
in dedicated low-intensity LED runs. The calibration 
is verified using singles collected continuously during 
regular data running. In either case, samples dominated 
by single photo-electron (SPE) hits are selected and fit 
after baseline subtraction using a function as shown in 
Eq.ffl 



S(x) 



- V tL 



l 



(x-nQif 

-_exp( — — ), 

r 2ncr, 



(4) 



which is a convolution of Poisson distributions with a 
Gaussian. Care is taken to ensure that SPE hits are dom- 
inant in the fit by requiring the PMT occupancy to be 
less than 0.13 and that the preceding analog baseline is 
stable. An example of the fitting is shown in Fig. [5] 

The fit values are stored in a database for use dur- 
ing data processing. We find a systematic difference 
in the average ADC count per SPE between AD1 and 




60 80 100 
ADC-preADC 



Figure 5: Fit to a single PMT SPE distribution using the function 
specified in equation [4] after subtracting the baseline. Dashed curves 
show the SPE and double PE components. The units are in raw ADC 
counts. 



AD2. We also observe a slight drift in the calibration 
as a function of time; however, all channels appear to 
drift together and the difference between AD1 and AD2 
is stable, as shown in Fig. [6] The distribution of cali- 
bration constants, ADC counts per photo-electron (p.e.), 
for AD1 and AD2 is shown in Fig. [7] for one of the cal- 
ibration periods. The channel to channel variation is 
within expectations. 



o DYB-AD1 -LED Calibration Gain 
o DYB-AD2-LED Calibration Gain 
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Figure 6: The stability of the peak in the SPE distributions averaged 
across each AD over a period of three months. 



4.2. PMT singles rate 

The PMT singles rates are studied using periodic trig- 
gers. The observed dark rates for AD1 and AD2 are 
relatively stable and slowly decreasing (by about 10%) 
over several months, as shown in Fig. [8] With the chan- 
nel thresholds set near 1/4 p.e., the average dark rates of 
all PMTs are about 10 kHz for AD1 and about 12 kHz 
for AD2. 
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Figure 7: Distribution of SPE peaks for AD1 and AD2 for data taken 
in late December. 
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Figure 9: Trigger efficiency as a function of reconstructed energy at 
the edge of the AD target volume (r = 120 cm, z = 135 cm). The top 
figure is for Nhit triggers while the bottom figure is for E-Sum trig- 
gers. Triangles represent efficiency measurements from 68 Ge source 
data. Circles result from LED scans. The colored lines show best fits 
based on error functions. The black lines indicate minimum recon- 
structed energies E m \ n of IBD positrons. 



Figure 8: Distribution of PMT dark rates for AD1 and AD2. 
4.3. Trigger 

The performance of the AD trigger system was 
probed by placing radioactive sources and LEDs at dif- 
ferent locations within the target volume. Positrons 
from a 68 Ge source were used as a proxy for IBD prompt 
events. The threshold of one of the two redundant trig- 
ger modes (Nhit and Esum) is lowered to test the trigger 
response of the other mode. In order to measure the 
trigger response to weak LED light, an external trigger, 
synchronized to the LED driver, is used. Both methods 
gave consistent results. 

Figure [9] shows the measured trigger efficiency as 
a function of reconstructed energy at the edge of the 
AD target volume. The energy-equivalent of the trigger 
thresholds, determined by fitting error functions against 
the data, was measured at E± ~ 0.37 MeV for both trig- 
ger modes. 



The energy reconstruction will be described in later 
sections and does not correct for the non-linear re- 
sponse of the scintillator. Thus, a positron annihilation 
is reconstructed at a mean energy of E m [ n ~ 0.9 MeV 
with an energy resolution cr E ~ 0.1 MeV. At energies 
£min - 3<X£, the trigger efficiency is still e = 1^002 
throughout the scintillating volume, thus ensuring neg- 
ligible trigger inefficiency. 

5. Comparison of rates and energy spectra 

5.7. Calibration sources 

The ADs are calibrated periodically using the au- 
tomated calibration units on the lid of each detector. 
The calibration program includes routine weekly cali- 
bration, as well as dedicated detector studies. During 
commissioning, the detector response is studied exten- 
sively with high statistics using all three sources in each 
of the three ACUs, including a fine-grain scan along 
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each vertical axis. During weekly calibration, the com- 
bined Am- 13 C/ 60 Co sources are lowered separately to 
five vertical locations for all three ACUs. Data are col- 
lected for five minutes at each position. This is followed 
by a 68 Ge source deployed to five vertical locations for 
ACU-A, each for five minutes of data taking. 

The energy spectra observed using the Am- 13 C/ 60 Co 
sources at the center of AD1 and AD2 are shown in 
Fig. [lOl Background spectra are taken from physics 
runs (without source deployment) and subtracted by 
normalizing the background run time to the calibration 
run time. Backgrounds can be reduced by applying a 
vertex cut around the source position. Besides the neu- 
tron peak at 8.05 MeV and the cobalt peak at 2.5 1 MeV, 
there is a single gamma shoulder at about 1.5 MeV and 
two peaks at lower energies. The spectra for AD1 and 
AD2 are in good agreement. The source activities mea- 
sured by fitting the spectra are similar for the two ADs. 



from the asymmetry of the spectrum. The small Gaus- 
sian component models the absorption of a gamma ray 
from the 60 Co decay by the source packaging material. 
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Figure 11: The cobalt spectrum fit using a Crystal Ball function plus 
a Gaussian. 
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Figure 10: The energy spectra observed using the Am- 13 C/ 60 Co 
sources at the center of AD1 and AD2 are overlaid. 

The cobalt peak is fit with a Crystal Ball function lEHl 
plus a Gaussian (see Fig.[TT]). The Crystal Ball function 
is given by 



f(x;a,n,x,cr) = N-{ 2of 



exp(-^#) for ^ > -a 



for — < -a 



(5) 



where: 



n 

B = --\a\ 
\a\ 



N is a normalization factor and a, n, x and cr are pa- 
rameters which are fit with the data. The Crystal Ball 
function is chosen to separate the Gaussian response of 
the detector from the low energy tail, thus avoiding bias 



The Am- 13 C source produces neutrons with an av- 
erage kinetic energy of 3.5 MeV. Although the rate is 
low, neutrons from the Am- 13 C source can be cleanly 
selected by requiring a coincidence between the prompt 
proton recoil signal and the delayed neutron signal. No 
background subtraction or vertex cut selection is re- 
quired. The prompt signal from the proton recoil is 
shown in Fig. [T2j The ADs have matching spectra, in- 
dicating that the energy response and quenching of the 
Gd-LS is the same in both. To better quantify the dif- 
ference in the observed spectra, Fig.[l2]plots a variable 
called Asymmetry defined as 



Asymmetry 



(Nadi +N Am )l2 



(6) 



where A^adi,2 is the bin content for AD1 or 2, respec- 
tively. 

When a neutron is captured by Gd, it has a 0.1848 
probability of being captured by 155 Gd releasing gamma 
rays with a total energy of 8.536 MeV, and a 0.8151 
probability of being captured by 157 Gd, releasing 7.937 
MeV. Contributions from other Gd isotopes can be ig- 
nored. For an event at the center of the detector, the 
neutron peak can be fit with two Gaussians with their 
relative area constrained by their neutron capture prob- 
abilities and the relative widths constrained by 1 / V# 
energy resolution, as shown in Fig. [I3j The free pa- 
rameter for the peak is taken as the average of the two 
peaks weighted by their neutron capture probabilities. 
For an event close to the edge of the IAV, two Crys- 
tal Ball functions are used for fitting, with similar con- 
straints as above. 
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Figure 12: The prompt signal from the proton recoil of Am- C neu- 
tron. 
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Figure 14: The neutron capture time on Gd from the Am- 13 C source 
at the detector center. 




Figure 13: The neutron spectrum of the Am- 13 C source at the detector 
center, fitted with two constrained Gaussians. 



The neutron capture time on Gd with the Am- 13 C 
source at the detector center is shown in Fig. [14] Since 
the capture time is directly related to the Gd concentra- 
tion, Fig. [14] provides a measure of the identicalness of 
the Gd concentration in AD1 and AD2. The measured 
capture times are 28.70 ± 0.15 [is and 28.60 ± 0.15 fis 
for AD1 and AD2, respectively. 

Fig. [I5]displays the observed spectrum from the 68 Ge 
source deployed at the center of the detector. 60 Co con- 
tamination in the source is responsible for the peak at 
2.5 MeV and the tail, which is the shoulder in Fig. [TT] 
extending to lower energy. While undesirable, the im- 
pact of the 60 Co contamination on the fitting of the 68 Ge 
peak is negligible. The 68 Ge activities and spectra for 
the two ADs are similar while the contaminations are 
different. The peak is fit with two Gaussians, one de- 



scribing the full energy deposit in the detector and the 
other describing the energy loss in the source enclosure, 
as shown in Fig. [16] 
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Figure 15: The energy spectrum of the Ge source. 



5.2. AD energy response 

The AD energy is reconstructed based on the total 
charge of an event, which is the charge sum of all 192 
8 -inch PMTs. The visible energy E Y { S is defined as the 
total charge divided by an energy calibration constant. 



Qtot/C 



(7) 



We use two independent energy reconstructions, as de- 
scribed in Sec. 15.31 One approach is to determine the 
calibration constant by constraining the energy peak of 
the 60 Co at the detector center to 2.506 MeV. This en- 
ergy calibration is done weekly, at the same frequency 
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Figure 16: The 68 Ge spectrum fitted with two Gaussians. 



as the PMT gain calibration. A second approach using 
spallation neutrons captured on Gd is also used to de- 
rive and monitor the energy calibration constant. For 
this method, the calibration constants are updated more 
frequently during regular physics data runs. 

For the 60 Co calibration, the energy constant is about 
162 p.e./MeV for AD1 and about 163 p.e./MeV for 
AD2. The time dependence is shown in Fig. [TTJ For 
each AD, the values vary within a narrow band with a 
width of 0.4%. The time dependence is also monitored 
using the visible energy from the 60 Co source at five ver- 
tical locations for all three ACUs, and the 68 Ge source 
from ACU-A deployed at five vertical locations. The 
time dependence of the reconstructed energy is largely 
corrected by the weekly calibration. The systematic 
uncertainty is determined by taking an average of the 
differences between consecutive calibrations (p.e. per 
MeV), which is 0.2%. 
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Figure 17: The time stability of the energy calibration constants. The 
drift relative to the first data point is shown in the bottom panel. 

For a perfect detector with reflective panels at the 
top and bottom, the energy response should be uniform 



along the vertical axis. A z-scan along the central axis 
using the 60 Co source is shown in Fig.|T8j A z-position 
of -150 cm corresponds to the bottom of the 3-m AV 
while +150 cm is the top. The detector response shows 
good uniformity along z direction; however, two fea- 
tures exist. One is a decrease in response at the top and 
bottom due to absorption in the reinforcing ribs of the 
AVs and the non-ideal reflectivity of the reflectors. The 
other is an asymmetry introduced by the reflector loca- 
tions. The PMT array is centered on the bulk of liquid 
scintillator; however, the 3-m and 4-m AVs have cone- 
shaped tops. Therefore, the gap between the top ring 
of PMTs and the top reflector is larger than the gap be- 
tween the bottom ring of PMTs and the bottom reflector 
by 77 mm. This is a geometric effect that is well re- 
produced in Monte Carlo simulations. If we center the 
PMT array between the reflectors, the simulated detec- 
tor response becomes symmetric along the z axis. The 
same non-uniformity is also verified by 68 Ge source and 
spallation neutrons, as shown in Fig. Q21 For the spal- 
lation neutrons, only those in a small cylinder around 
the central axis are selected for comparing with ACU-A 
data. 
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Figure 18: The AD energy response as a function of z for Monte Carlo 
data, spallation neutron data, and ACU-A sources. 



The z-scan with the 60 Co sources from all three ACUs 
is shown in Fig. [19] The energy response for ACU- 
B and ACU-C exhibit similar features as ACU-A de- 
scribed above, but with more visible energy. This is 
also a geometric effect caused by an increase in aver- 
age acceptance as the source approaches the PMTs. It 
is well reproduced in Monte Carlo, and is linear in r 2 , 
where r is the radial position of the source. The dif- 
ferences between the ADs are small for all three ACU 
scans. Analogous to Eq. [6) we define the asymmetry 
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between the fitted energy peaks {ji) 



Asymmetry = 



G^ADl +y^AD2)/2 



(8) 



The Asymmetry values lie within a narrow band with a 
width of 0.3% for ACU-A, 0.4% for ACU-B, and 2% 
for ACU-C. The ACU-C scan is in the gamma catcher, 
and may suffer more from edge effects, such as acrylic 
properties and PMT-to-PMT efficiency variations. 
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Figure 19: Visible energy in z with respect to the detector center for 
60 Co sources from each ACU is shown at top. The corresponding 
Asymmetry parameter is shown for ACU-A (Asy. A), ACU-B (Asy. 
B) and ACU-C (Asy. C) in the bottom panels. 

Detector response (observed energy versus deposited 
energy) is not linear in energy, and the effect differs with 
particle species. This nonlinearity is due to the quench- 
ing effects of the liquid scintillator, Cherenkov light 
emission and subsequent absorption and re-emission in 
the liquid scintillator, PMT dark noise, PMT or elec- 
tronic non-linearity, etc. However, the effects are the 
same in each AD as shown in Fig. [20l In this figure, 
the three radioactive sources are deployed at the detec- 
tor center, and all alphas, spallation neutrons, and IBD 
neutrons are distributed uniformly throughout the detec- 
tor. The alphas result from daughter decays from resid- 
ual contamination of 238 U, 232 Th and 227 Ac. Clean sam- 
ples of alphas from 214 Po, 212 Po and 215 Po decays are 
selected using the ft - a or a - a time correlations. The 
energy is vertex-corrected, as described in Sec. 15.31 

5.3. Energy reconstruction 

Given the importance of minimizing the energy scale 
uncertainties in the Daya Bay analysis, two indepen- 
dent energy reconstructions are adopted to cross check 
each other. One is based on source calibration data and 
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Figure 20: Asymmetry of energy nonlinearity for two ADs. See text 
for details. 



the other uses spallation neutrons. Both apply a ver- 
tex correction to the visible energy. Since minimizing 
the event selection criteria is important to maintaining a 
small systematic error, vertex reconstruction will not be 
used to select IBD events. However, vertex reconstruc- 
tion is necessary to correct the position dependent en- 
ergy non-uniformity. It is also valuable for many other 
studies. 

5.3.1. Energy reconstruction based on sources 

The vertex reconstruction is based on center-of- 
charge (COC), defined as the charge- weighted-mean of 
the coordinates of all PMTs. The mapping from COC 
to vertex is done by analytic corrections. The correction 
formula is inspired from Monte Carlo simulations and 
is almost linear in COC but with small corrections. The 
reconstructed vertex is 



c\ x Rqoc ~c 2 x R coc , 



(9) 



z = (Zcoc - c 3 x Z coc ) x (c 4 - c 5 x Rcoc) , (10) 

where c\ - c$ are coefficients that are found by calibrat- 
ing the COC to the true locations of the cobalt sources 
with data. Rcoc and Zcoc are the COC coordinates. 

To correct the vertex-dependent energy non- 
uniformity, cobalt source z-scan data is used to model 
the corrections in terms of the radius r and z. The 
reconstructed energy is 



7 Co 



C n Co x£ vls /(C r Co (r)xC 7 Co (z)), 



(11) 



where E Y { S is defined in section 5.2, Cjr°(r) is a linear 
function of r 2 , Cf°(z) is a 3rd-order polynomial, and 
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C„° is a constant determined by the ratio of the neu- 
tron peak of the Am- 13 C source to the cobalt peak, both 
at the detector center. C„° rescales the corrected energy 
from the cobalt energy scale to the neutron energy scale. 
Both AD1 and AD2 data are used to extract the coeffi- 
cients in Eq. l9lfTT1 All AD events will be reconstructed 
to the neutron energy scale first. 

As displayed in Fig. [21] after the vertex-dependent 
correction, the non-uniformity of visible energy shown 
in Fig. [19] has been minimized. All data points for each 
ACU in both ADs now lie within a band from 0.98 to 
1.02, relative to the cobalt energy response at the detec- 
tor center. Since we did not apply detector-dependent 
corrections, and the vertex-dependent correction fac- 
tor C^°(r) x Cf°(z) is at most 10% within the target 
and gamma catcher volume, the reconstructed energy 
has nearly the same relative non-uniformity between the 
ADs as the visible energy. 
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Figure 21: Reconstructed energy in z with respect to the detector cen- 
ter for 60 Co sources from each ACU, after the vertex-dependent cor- 
rection. 



5.3.2. Energy reconstruction based on spallation neu- 
trons 

More than 8,000 spallation neutrons per detector per 
day are observed at the near site. They can be cleanly 
selected by searching in a time window from 20 to 200 
jjs after a muon passes through an AD. The lower limit 
is required to ensure that the baseline of the electronics 
has recovered after a large muon signal. Capture on both 
Gd and H can be seen in the spallation neutron spectrum 
in Fig.[22l The small Asymmetry shows that the spectra 
for the two ADs are almost the same. 

The vertex position is reconstructed by comparing the 
number of photo-electrons observed in each PMT with 
templates produced through MC simulation. Templates 



are produced for 20 bins in the r direction, 20 bins in the 
z direction, and 24 bins in the direction. For each grid 
point where the templates are produced, we compute the 
X 2 values defined as: 



X 



PMTs 

2 



W bs ,A^ xp (r,z,0)) 



-2 In 



P(N obs ,N obs ) 



(12) 



where P(n,fi) = id n e~^ /(n\) is the probability of find- 
ing n photo-electrons when the mean value is \i assum- 
ing Poisson statistics, N° hs is the observed number of 
photo-electrons and A^ xp (r,z, 0) is the expected num- 
ber of photo-electrons from the templates. The recon- 
structed vertex position is calculated by interpolating 
the^f 2 distribution. 

To correct the vertex-dependent energy non- 
uniformity, spallation neutron data is used to model the 
corrections in terms of radius r and z. The reconstructed 
energy is given by: 



7 SN 



C s r N (r)xC s z N (z)x Qt0 JC s e N , 



(13) 



where C s r N (r) and C% (z) are 3rd-order polynomials de- 
rived from the observed non-uniformity of the spalla- 
tion neutron energy response in the AD, gtot is the ob- 
served total charge, and Cf N is a charge to energy con- 
version coefficient determined by constraining the spal- 
lation neutron peak to 8.047 MeV. 
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Figure 22: Spallation neutron energy spectrum and Asymmetry be- 
tween ADs. 

The charge to energy conversion coefficient C s e N 
is approximately 168 p.e./MeV for AD1 and 169 
p.e./MeV for AD2. The variation of this calibration 
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constant with time is shown in Fig. [23j Although data 
from special runs are included in the plot (e.g. PMT sys- 
tematic studies) data from these special runs are not in- 
cluded in the side-by- side comparison. Excluding these 
runs, the absolute variation is less than 1%. The time 
dependence of the energy reconstruction is estimated by 
monitoring the 60 Co calibration data. The reconstructed 
energy of the 60 Co sources at the detector center varies 
0.3% for both ADs, while the rms of the Asymmetry 
drift over time is only 0.15%. The correction to the 
non-uniformity is also checked with the 60 Co data, as 
shown in Fig. [24l After vertex correction, the energy 
non-uniformity is almost the same as in Fig. [211 
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Figure 23: Variation in the spallation neutron energy conversion coef- 
ficient with time is shown at top. The associated Asymmetry param- 
eter is shown at bottom. The dashed vertical lines bound a period of 
dedicated systematic studies. 
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Figure 24: Variation in the 60 Co energy peak with source location 
after spallation neutron based vertex correction. 

The energy reconstruction methods are consistent 



with each other, and each offers distinct advantages. 
The source data method has the advantages of extremely 
low background, high statistics, and identicalness at 
near and far sites. Furthermore, the source calibration 
data alone is sufficient to reconstruct the spatial uni- 
formly distributed IBD neutron peak and the spallation 
neutron peak at their true energy, as we will see in later 
sections. This implies that the detector response is well 
understood. 

Spallation neutrons have a spatial distribution nearly 
as uniform as the IBD neutrons. Therefore, this energy 
calibration method equally samples the entire target vol- 
ume and corrects the spatial non-uniformity of IBD neu- 
trons more accurately. Furthermore, spallation neutron 
events are extracted from all regular physics runs, com- 
pared to the weekly source calibrations. 

The resolution of the reconstructed energy is shown 
in Fig.O The resolution curve is fitted phenomenolog- 
ically to be (7.5/ V^(MeV) + 0.9)%. 
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Figure 25: Resolution of reconstructed energy. 



5.4. AD triggers 

AD data is primarily collected through an OR of the 
multiplicity and energy sum triggers with an energy 
threshold around 0.4 MeV (see Fig. 9). The typical 
trigger rate is about 280 Hz per AD. About 5% of these 
are so-called flasher events, instrumental background 
events apparently resulting from an electronic discharge 
in the PMT assembly. The observed energy of flasher 
events ranges from threshold to 100 MeV. About 5% 
of the PMTs have been identified as sources of flasher 
events. The flasher events have a geometric charge pat- 
tern that is distinct from other physical events. As such, 
they can be easily identified and removed. Fig. [26l shows 



15 



the discrimination of flasher events for the delayed sig- 
nal of IBD candidates, where the discrimination vari- 
able FID > indicates a flasher. For this demonstration, 
flasher events have been kept in the IBD selection. For 
the IBD analysis as well as most other analyses, the re- 
jection of flasher events will be done at the beginning 
of the data reduction. Most of the events around FID 
= are actually pile-up events that confuse the discrim- 
ination algorithm. The distributions for AD1 and AD2 
overlay well for the selected IBD events (FID < 0), but 
there is some variation between the distributions when 
FID > 0. Monte Carlo studies indicate that the FID is 
extremely efficient at identifying flasher events. The in- 
efficiency for selection and contamination of IBD selec- 
tion are evaluated to be 0.02% and < 10" 4 respectively, 
for both AD1 and AD2. The relative uncertainties are 
estimated to be 0.01%. Special runs were conducted to 
identify PMTs that exhibit flashing; however, due to the 
high efficiency of the FID, 100% of the AD PMTs are 
operational, including the PMTs identified as flashers. 




Figure 27: Spectrum comparison of all AD triggers for AD1 and AD2 
after flasher rejection. The lower plot measures the Asymmetry be- 
tween AD1 and AD2, which is largely due to nuisance triggers imme- 
diately following a muon event. 




FID of delayed signal 

Figure 26: Flasher event discrimination for the delayed signal of IBD 
candidates. Events with a PMT discharge (flasher) have a FID > 
while real IBD events have a FID < 0. 

The spectra of AD triggers for both detectors are 
shown in Fig. [27] after removing flasher events, along 
with the associated Asymmetry parameter. The knee at 
~10 3 MeV corresponds to the maximum path-length by 
a minimum ionizing muon through the 4-m diameter by 
4-m high active detector volume. Higher energy events 
are showers induced by muons. The greatest difference 
between the spectra, about -15%, is created in large 
part by triggers that occur shortly after a muon passes 
through a detector. These nuisance triggers may last up 
to 10 yus following a high-energy muon, and result from 
PMT after-pulsing and ringing, and signal overshoot in 
the readout chain i22ll . 

The performance of the Day a Bay muon system will 
be described elsewhere. The event rates for a multi- 
plicity of twelve or more PMTs in the IWS and OWS 



are shown in Fig. [28J There are 121 PMTs installed in 
the IWS and 167 PMTs in the OWS. Occasionally, the 
muon rates are disturbed by electronics noise. The ob- 
served rate of muons in each AD is 21 Hz, assuming that 
any AD trigger with an energy >20 MeV results from a 
muon. For such an AD muon, the detection efficiency 
of the IWS is 99.7%, and that of the OWS is -97%. 
We believe that these are lower limits since high energy 
AD events may be due to physical processes other than 
muons (e.g. from high energy neutrons produced in the 
surrounding rock). Such events have been observed in 
Monte Carlo simulations. The apparent inefficiency of 
the OWS is largely due to muon decay in the IWS or an 
AD. 
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Figure 28: Muon rate in the IWS and OWS with a PMT multiplicity 
>12. 

Figure[29]displays the AD spectra after excluding AD 
triggers in the time interval (-2, 200) jjs with respect 
to an IWS or OWS trigger. The -2 yus cut is used to 
avoid time alignment issues between detectors. Besides 
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the removal of AD events with muons or muon daugh- 
ter products, nuisance triggers that closely follow AD 
muons are also rejected. Good agreement is observed in 
the spectra between the ADs. 
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Figure 29: Spectrum comparison of AD triggers after muon veto cuts 
have been applied. The associated Asymmetry distribution is also 
shown. 



5.5. IBD rates and energy spectra 

To select IBD events, the following criteria are ap- 
plied: 

1 . the flasher rejection cut described above; 

2. positron energy: 0.7 MeV < £ pr0 mpt < 12.0 MeV; 

3. neutron energy: 6.0 MeV < delayed < 12.0 MeV; 

4. time coincidence: 1 yus < Ar prom p t 4 e i aye d < 200 yus; 

5. a muon veto that rejects a prompt-delayed pair 
if the delayed signal is within 600 /is after an 
IWS or OWS trigger; however, if a muon deposits 
> 20 MeV in an AD the exclusion window is ex- 
tended to 1000 yus due to the increased probability 
of multiple neutrons. If a muon deposits > 2.5 GeV 
in an AD, the exclusion window extends to 1 .0 sec- 
ond to reject long-lived cosmogenic backgrounds. 

6. a multiplicity cut that requires no other > 0.7 MeV 
trigger 200 yus before the prompt signal and 200 yus 
after the delayed signal; 

A scatter plot of prompt-delayed energy pairs is 
shown in Fig.[30| In the plot, IBD candidates with neu- 
tron capture on hydrogen can also be identified although 
the low energy part is obscured by accidental coinci- 
dence background. 
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Figure 30: Scatter plot of the prompt-delayed signals. Pairs inside the 
dotted line box are taken as IBD candidates. 

Two detector-dependent inefficiencies, the muon veto 
dead time and the multiplicity cut, have been corrected 
in the following studies on the IBD rate and spectrum 
comparison. The energy selection, time coincidence se- 
lection, as well as other absolute IBD efficiencies such 
as the fraction of neutron capture on Gd are not cor- 
rected. They are identical for ADs from our analysis 
thus not relevant to the side-by-side comparison of ADs. 
However, the uncertainties are analyzed. 

All triggers within 2 /is before an IWS or OWS muon 
are also vetoed to avoid time alignment issues between 
detectors. IWS and OWS veto time intervals are care- 
fully merged with the veto time interval of each AD to 
avoid double counting. All livetime segments between 
muons are precisely calculated. The total veto dead time 
is measured to be 17.68% for AD1 and 18.06% for AD2 
where the difference comes from the AD response to 
high-charge events and resulting muon rejection.. 

The multiplicity cut efficiency consists of three parts. 

1. No other > 0.7 MeV trigger (singles) in a 200 jis 
window before the prompt signal. The efficiency is 
61 = 1 - R s • 200 /is where R s is the rate of singles, 
which is ~60 Hz in both ADs. 

2. No singles between the prompt and delayed sig- 
nals. The efficiency is 62 = 1 - R S T C , where T c is 
the average neutron capture time. 

3. No singles in a 200 /is window after the delayed 
signal. The efficiency is 63 = l-R s - 200 jis(l - 100 
jds/T v ) for a single livetime segment T v longer than 
200 yus, or 63 = 1 - R s T v /2 for T v shorter than 200 
yus. 63 is different from e\ because a muon may 
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occur and veto the following singles, thus reducing 
the rejection probability. 

The multiplicity cut efficiency is evaluated using Pois- 
son statistics and the probability density function of the 
neutron capture time. The difference from the above 
approximation is < 0.015%. The inefficiency from the 
multiplicity cut is -2.5%, depending on the singles rate. 
Combining the multiplicity cut and the muon veto, the 
total efficiency for these two cuts are 80.10% for AD1 
and 79.76% for AD2. 

Another analysis has been done with a different mul- 
tiplicity cut that requires: 

1 . Only one prompt candidate within 200 yus before a 
neutron candidate; 

2. No other prompt candidate within 400 yus before 
the neutron candidate; 

3. No other neutron candidate within 200 yus after the 
first neutron candidate; 

4. At the same time, applying an additional muon 
veto 200 yus before each IWS or OWS muon. 

The last cut is required to avoid a high energy (> 6 MeV) 
prompt signal of an IBD being mistaken as a neutron 
candidate and forming an accidental background with 
a preceding background gamma while the real neutron 
signal is vetoed by a muon. This multiplicity cut has 
a fixed time window and thus is decoupled from the 
neutron capture time and the muon veto, with the cost 
of dropping ~4% more IBD events. The selection ef- 
ficiency is exp(-R s • 400 /us) exp(-7?j • 200 /us) where 
Rd is the rate of neutron candidate. The two multiplic- 
ity cuts have been compared and verify the efficiency 
calculations. 

For the purpose of this side-by-side comparison, we 
have evaluated the three largest backgrounds in the IBD 
sample: accidentals, 9 Li/ 8 He decays, and fast neutrons. 
Other correlated backgrounds such as (a,n) reactions or 
the ACU sources have been evaluated and found to be 
negligible. 

The accidental background rate was determined by 
separately counting the singles rate for both £ + -like and 
neutron-like signals. We get 10.20 ± 0.05/day for AD1 
and 10.10 ± 0.05/day for AD2, corresponding to -1.7% 
Background- Signal ratio (B/S) and a 0.01% uncertainty 
in the analyzed data set. 

The correlated background from the /3-n cascade of 
9 Li/ 8 He decays is evaluated by fitting the time between 
the last muon and prompt IBD candidate l23ll with the 
known decay times for these isotopes. The 9 Li/ 8 He de- 
tection rate is found to be 4.2 ± 1.2/day. The 9 Li/ 8 He 
detection rates for muons that deposit < 2 GeV in an 



AD is consistent with zero within statistical error. After 
applying the 1 -second shower muon veto, the 9 Li/ 8 He 
background in the IBD sample is estimated to be < 0.3% 
(1 cr). 

The fast neutron backgrounds are estimated in two 
ways. 

1 . The upper limit of the prompt energy cut is relaxed 
to extend the prompt energy spectrum to high en- 
ergy. A flat distribution is observed for energies 
> 12 MeV. Assuming the energy spectrum for fast 
neutrons is flat through the relevant neutrino en- 
ergy region, we estimate the relative rate of the fast 
neutron background within the IBD sample to be 
0.2%. 

2. For muons tagged by the IWS and OWS, we get the 
rate and energy spectrum of the fast neutron back- 
grounds as a function of muon track length in the 
water pool. A fast neutron can leak into the IBD 
sample if its parent muon is not detected due to in- 
creased inefficiency for detecting muons of short 
track length in the water pool. Estimates of the fast 
neutron backgrounds produced by muons passing 
through nearby rock rely on simulations. The B/S 
is measured to be 0. 15 ± 0.05%, consistent with the 
first method. 

For the side-by-side comparison of AD1 and AD2 in 
the same water pool, the background contents are very 
similar thus not critical for the relative comparison. 

The prompt energy spectra of selected IBD events 
are shown in Fig. [3J] for AD1 and AD2. The back- 
grounds, dominated by accidental coincidence, are sub- 
tracted statistically in the plot. The ratio of the total IBD 
rates in AD1 and AD2 is 0.987 ± 0.008 (stat). 

The neutron energy peak is a critical check on the 
energy scale calibration and related uncertainties. The 
neutron energy distribution of the selected IBD sam- 
ples is shown in Fig. [32] When fitted with double 
Crystal Ball functions described in Sec l5.1[ the neutron 
peak agrees very well with the expected value of 8.047 
MeV and the expected resolution. The peak for AD1 is 
~0.3% higher than that for AD2, and will be discussed 
in Sec. 16.21 

The capture time spectrum has four components. The 
neutrons from the IBD reactions thermalize in ~8 /us 
forming the rise in the spectrum. Thereafter, the cap- 
ture on Gd dominates and forms the exponential com- 
ponent. At times larger than 100 /us, the spill-in/spill- 
out effect^] at the edge of the target vessel play an im- 
portant role. Neutrons produced in the gamma-catcher 

1 'Spill-in' ('spill-out') occurs when an IBD neutron produced out- 
side (inside) the Gd-LS volume is detected in (outside) the Gd-LS 
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Figure 3 1 : Energy spectrum for the prompt signal of IBD events in 
each AD along with the ratio in each energy bin. The dashed line 
shows the ratio R of the total rates in ADl and AD2. 

survive for longer periods of time, allowing them an op- 
portunity to spill-in to the target volume before being 
captured on Gd. Finally, there is a small flat component 
from accidental coincidences. The neutron capture time 
on Gd for the whole IBD sample is shown in Fig. [33] 
The expected accidental background component is also 
shown. To compare with the Am- 13 C results, the sam- 
ple is purified by applying a vertex cut (r < 1.25 m and 
\z\ < 1.25 m) on the delayed signal events to reject spill- 
in candidates, and a cut on prompt energy (£ pr0 mpt > 3 
MeV) to further reject the accidental backgrounds. The 
neutron capture times for the cleaned IBD samples are 
28.2+0.3 yus for ADl and 28.6+0.3 //s for AD2, in good 
agreement with the Am- 13 C results. 

The IBD rates are shown in Fig. [34] as a function 
of time. The D2 core was shut down on October 26, 
2011 and returned to service on December 9. Nor- 
mally, reactor cores are ramped-up over 1-2 weeks be- 
fore reaching full power. At full power, the D2 core 
contributes around 40% of the neutrinos observed by 
ADl and AD2. The ramping can clearly be seen in the 
observed IBD rate. On November 9, the L3 core came 
back online after refueling. On December 12, the L2 
core was shutdown. The L2 and L3 cores contribute 
-8% and ~3% of the neutrinos observed by ADl and 
AD2, respectively. During the three-months of data tak- 
ing, there were two periods (shown as vertical shaded 
areas) amounting to 203.2 hours that were dedicated to 



volume. These phenomena effectively alter the target mass for an- 
tineutrinos. 



Figure 32: IBD neutron energy distribution for each AD along with 
the associated Asymmetry distribution. 

test runs (e.g. flasher and PMT testing). 

6. Detection uncertainties 

6.1. Target mass 

The target mass is defined as the mass of Gd-LS in- 
side the 3-m inner acrylic vessel (IAV), up to the top of 
its conical lid. Two sets of bellows extend above the 
IAV, connecting the IAV to the central overflow tank 
and a fluid monitoring camera. The target mass is de- 
termined from the total mass of Gd-LS added to the de- 
tector during filling, minus the mass of the fluid in the 
overflow tanks and the two bellows. The fluid height in 
the overflow tank is monitored by a multiply-redundant 
system of sensors. The mass of Gd-LS in the bellows 
is computed from the bellows volume and the Gd-LS 
density. 

The total number of free protons (hydrogen nuclei) is 
simply given by the target mass multiplied by the num- 
ber of hydrogen nuclei per kg of Gd-LS. The Gd-LS 
mass density is described by the formula 

p t =po/a+J3(T-T )), (14) 

where p = 0.8590 kg/L at T = 19 C and p = 
9.02 + 0.16 x 10" 4 . Based on a purely analytical calcu- 
lation incorporating a realistic distribution of LAB hy- 
drocarbon chains, the proton mass fraction in Gd-LS is 
determined to be 11.77%. A combustion-based mea- 
surement of the hydrogen fraction in Gd-LS is used in 
this analysis and gives 12.01+0.42% and 11.97 + 0.47% 
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Figure 33: Neutron capture time for the selected IBD sample along 
with the associated Asymmetry distribution. 
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Figure 34: IBD rate in ADl and AD2 as a function of time. 



in ADs 1&2, respectively, including an instrumental un- 
certainty of 0.3%. Since the measurements and calcula- 
tion are consistent, and the ADs are filled equally from 
five Gd-LS storage tanks, we consider the relative un- 
certainty on the free proton mass fraction between de- 
tectors negligible. For the absolute mass fraction uncer- 
tainty, we take the largest uncertainty, 0.47%. 

There are no ribs or other structures inside the 3m 
vessel to trap bubbles. The IAV is a vertical right cylin- 
der with inner diameter 3.10 m and inner height 3.07 
m. The cylinder is topped by a shallow cone with di- 
ameter 3.04 m and acute 3 degree edge angle. The total 
volume is 23.364 m 3 . Multiplying by a density of 860 
kg/m 3 gives a nominal target mass of 20,093 kg. Since 
the blinding is 120 kg, a nominal target mass of 20,000 
kg is assumed for simplicity. 

The total mass of Gd-LS inside the AD is determined 
during filling. This quantity is subsequently blinded at 
the level of 0.6% (120 kg). The mass of Gd-LS is de- 
termined primarily from load cells under the ISO tank, 



with a backup measurement from a Coriolis flow me- 
ter. Multiple calibrations during ISO tank filling indi- 
cate that the four load cells are linear. Long-term tests 
with the load cells show that the best precision is ob- 
tained from a 15-minute average of load cell readings, 
and that the load cell readings may drift by up to 3 kg 
over several hours. Corrections are made to account 
for the weight of dry-nitrogen that enters the ISO tank 
during filling. The mass as determined by the Coriolis 
meter is consistent within the accuracy of the meter for 
both ADs. Overall, the relative uncertainty on the to- 
tal mass is dominated by the irreducible load cell drift. 
The combined uncertainty on the total mass is 3.0 kg or 
0.015%. 

The overflow tanks at the top of the AD accommo- 
date fluid expansion and contraction from temperature 
changes. The mass in the overflow tanks is computed 
from the liquid levels, the tank geometry, and the liquid 
density. The liquid level is determined from a redundant 
set of lid sensors. The most accurate sensor uses ultra- 
sonic sound to measure the distance from the sensor to 
the fluid with 1 mm accuracy. Cross-checks on this sen- 
sor are provided by a capacitance sensor, and by camera 
observation of the liquid levels in the off-center calibra- 
tion port. Although the ultrasonic sensors have a res- 
olution of 0.1 mm, calibration reveals deviations from 
linearity at the 1-mm level. A 1 mm uncertainty in liq- 
uid height corresponds to a 1 .33 L, or 1 . 14 kg (0.0057%) 
uncertainty on the target mass. The conversion from liq- 
uid height to liquid volume is calculated from a survey 
of the overflow tank geometry, and was cross-checked 
by filling one of the tanks with deionized water in 2L 
increments. The largest difference between calculated 
and measured liquid volume is 1.5 L. Multiplying by a 
fluid density of 0.86 kg/L gives an uncertainty of 1.3 kg. 
Surveys of all overflow tanks show good similarity, so 
we use a general function to compute the overflow mass 
of all tanks with a maximum deviation from any physi- 
cal tank of 0.22 kg. Combining these two uncertainties 
gives a total uncertainty on the calculation of fluid mass 
from fluid height of 1.32 kg or 0.0066%. A tilted over- 
flow tank would cause an error in fluid height as mea- 
sured by the offcenter ultrasonic sensor. The levelness 
of the ADs has been measured, and the uncertainty from 
the tilt is less than 1.37 kg, or 0.0068%. 

Two identical bellows connect the 3m inner acrylic 
vessel to the SSV lid. The bellows are corrugated, 
with cylindrical cuffs on each end. The bellows do 
not compress during AD assembly, but the cuffs are 
somewhat free to slide in their housing. Assuming rela- 
tively large uncertainties on the bellows dimensions, 5% 
cross-section and 8% (6 cm) length, the bellows volume 
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is 4.30 ± 0.42 L or 3.70 ± 0.36 kg. The bellows are 
attached to the IAV lid at two stubs. The total fluid vol- 
ume in the two stubs is 5.78 L. The volume uncertainty 
is assumed to be equivalent to the bellows uncertainty, 
0.4 L. Overall, the bellows and stub uncertainty is 0.58 
L or 0.5 kg (0.0025%). The uncertainties on the number 
of free protons are summarized in Table [3j 



Table 3: Uncertainties on the number of target protons. Relative and 
absolute uncertainties are given for a single AD. Fractional uncertain- 
ties are computed assuming a nominal 20 t target mass. 



Quantity 


Relative 


Absolute 


Free protons/kg 


neg. 


0.47% 


Density (kg/L) 


neg. 


0.0002% 


Total mass 


0.015% 


0.015% 


Overflow tank geometry 


0.0066% 


0.0066% 


Overflow sensor calibration 


0.0057% 


0.0057% 


Overflow tank tilt 


0.0068% 


0.0068% 


Bellows Capacity 


0.0025% 


0.0025% 


Target mass 


0.019% 


0.019% 


Free protons 


0.019% 


0.470% 



6.2. Energy scale uncertainty 

As shown in Sec. 5.1, the energy scale uncertainty 
from the time variation is estimated to be ±0.2% by 
averaging the time variation of all ACUs. The non- 
uniformity differences between two ADs, when aver- 
aging ACU-A and ACU-B, lead to an uncertainty of 
±0.2%. The non-linearity from the cobalt energy scale 
to the neutron energy scale should be the same for the 
two ADs. We find that the difference is 0.1%. From 
the comparison of calibration data, we estimate that the 
relative energy scale uncertainty for neutrons uniformly 
distributed in the target volume is -0.3%. 

Besides the calibration data, the energy scale uncer- 
tainty related to the non-uniformity can be evaluated 
from the IBD neutron and spallation neutron distribu- 
tions. Fig. [35] displays the Asymmetry of energy peaks 
between the two ADs for both IBD and spallation neu- 
trons captured on Gd. The target volume is divided into 
twenty bins (pixels), four along R 2 and five along the z 
direction. The indexing of the bins is shown in the fig- 
ure. Each bin has equal volume thus equal numbers of 
neutrons, assuming the spatial distribution is uniform. 
The actual relative bin content could be slightly differ- 
ent due to the spill-in/spill-out effects, slowing down of 
spallation neutrons, vertex reconstruction and event se- 
lections. Events reconstructed outside the target volume 
have been included into the closest bin. The Asymme- 
try for IBD neutrons ranges from -0.1% to 1.0%, with a 



fitting error of 0.1%. Spallation neutrons have the same 
structure as IBD neutrons, verifying the non-uniformity 
differences between two ADs. The mean values of the 
Asymmetry in the twenty bins are 0.3% for both IBD 
neutrons and spallation neutrons. By the definition of 
the Asymmetry, it means that the source-based calibra- 
tion results in a 0.3% higher energy scale for AD1 than 
AD2, agreeing with the aggregate energy peak differ- 
ence shown in Fig. [32] and Fig. [36] The difference ap- 
pears in the upper half of the ADs. For the same rea- 
son, the energy Asymmetry of any uniformly distributed 
sources, such as alphas in Fig. [20] will be 0.3% higher 
than the Asymmetry of sources at the detector center. 
The RMS value of the bins is 0.27% for IBD neutrons 
and 0.21% for spallation neutrons, which will be taken 
as the relative energy scale uncertainties from the non- 
uniformity. 
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Figure 35: Asymmetry distribution in the peak of the energy distribu- 
tion shown bin by bin for spallation and IBD neutrons. 

The direct comparison of the IBD neutron energy has 
been shown in Fig. [32] The difference between the two 
ADs is 0.3% with AD1 being higher. For spallation neu- 
trons, the same comparison is shown in Fig. [36j The 
difference between the two ADs is also 0.3% with AD1 
being higher, but systematically higher than IBD neu- 
trons also by 0.3%. Combining the time variation, the 
non-uniformity, and the non-linearity uncertainties, we 
estimate that the relative energy scale uncertainty for 
the 8.047 MeV neutron peak for this pair of ADs is 
0.4%. This calibration-related uncertainty is reflected in 
the 0.3% relative energy scale difference observed after 
source-based calibration. 

The efficiency of the 6 MeV cut for the neutron se- 
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Figure 36: Spallation neutron energy distribution for each AD along 
with the associated Asymmetry. 

lection (Sec. 5. 5) suffers from edge effects in the ADs. 
The neutrons at the detector center are not affected by 
the 6 MeV cut selection since the energy is fully con- 
tained, as shown in Fig. \13\ When approaching the 
detector edge, the non-uniformity Asymmetry between 
the ADs increases. We use a high statistics IBD sim- 
ulation to generate a weight map corresponding to the 
non-uniformity map in Fig.[35l by counting the neutrons 
around 6 MeV. After re- weighting, the relative energy 
scale uncertainty for the 6 MeV cut is estimated to be 
0.45% for this pair of ADs. The relative efficiency un- 
certainty for the 6 MeV neutron energy cut, when the 
energy scale varies by 0.45%, is 0.1 1%. 

6.3. Summary of detector-related uncertainties 

Relative uncertainties from other analysis cuts ap- 
plied in the IBD selections are summarized below. 

1. The flasher rejection cut is rather safe and robust 
for the IBD selections as shown in Sec. 15.41 The 
uncertainty is estimated to be 0.01%. 

2. The AD trigger has almost 100% efficiency at 0.7 
MeV, which is 3<x away from the e + IBD threshold. 
The IBD reaction in acrylic may have prompt en- 
ergy lower than 0.7 MeV. Assuming a 2% energy 
scale uncertainty to account for the edge effect, we 
get 0.01% uncertainty for the £p r0 mpt selection effi- 
ciency by Monte Carlo simulation. 

3. In principle, the neutron capture time for the two 
ADs should be identical due to the specially de- 
signed Gd-LS mixing and filling procedure. From 



the fitted Am-C neutron capture time, IBD neu- 
tron, and spallation neutron capture time, the Gd 
concentration difference is estimated to be < 0.5%. 
By varying the Gd concentration in Monte Carlo, 
we find the uncertainty of the time coincidence cut 
to be < 0.01%. 

4. Of the three components of the multiplicity cut, e\ 
and 63 can be accurately calculated, up to the sta- 
tistical precision of singles. 62 relies on the aver- 
age neutron capture time. If the Gd content in the 
ADs is the same, the efficiency uncertainty is neg- 
ligible. By measuring the average neutron capture 
time, this uncertainty is estimated to be < 0.01%. 

5. From neutron capture time of Am- 13 C sources and 
IBD neutrons, we limit the relative fraction of neu- 
tron captures on Gd in the ADs to < 0.1%. 

6. Relative differences in acrylic vessel thickness and 
density, and liquid density between ADs can re- 
sult in relative differences in the number of n-Gd 
captures originating from outside the target region. 
Using MC simulation, the measured sub-percent 
relative differences in these quantities in AD1 and 
AD2 lead to an expected relative spill-in/out uncer- 
tainty of 0.02%. 

7. The determination of the livetime includes pre- 
cision calculations of the muon veto time and 
blocked triggers. The precision of the livetime is 
estimated to be < 0.01%. 

The detector-related relative uncertainties for AD1 
and AD2 are summarized in Table |4j The efficiencies 
of the multiplicity cut and capture time cuts, the relative 
precision on H/Gd ratio, and the relative uncertainties 
of spill-in/out effects are correlated; however, the H/Gd 
ratio uncertainty is dominant among the four. The total 
uncertainty is estimated to be 0.2%. 



Source of uncertainty 


Quantity 


Mass measurement relative precision 


0.02% 


Flasher cut 


0.01% 


Efficiency of neutron energy cut 


0.11% 


Efficiency of e + threshold cut 


0.01% 


Efficiency of multiplicity cut 


< 0.01% 


Efficiency of capture time cuts 


0.01% 


Relative precision on H/Gd ratio 


<0.1% 


Relative uncertainty of spill-in/out 


0.02% 


Livetime precision 


<0.01% 


Total detector-related uncertainty 


0.2% 



Table 4: Detector-related relative uncertainty of Daya Bay evaluated 
with the 1st pair of ADs in the Daya Bay Near Hall. 
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When the analyses were frozen, the baselines, the 
thermal power histories of the cores, and the target 
masses of the two ADs were unblinded. The expected 
ratio of the IBD rates in AD1 and AD2 is 0.981 as 
compared to the measured ratio 0.987 ± 0.008(stat) ± 
0.003 (sy st). The deviation from unity of this ratio is 
largely due to small differences in the baselines of the 
two ADs. The reactor flux differences can be ignored, 
and target mass difference contributes 0.15% to the de- 
viation. 

7. Conclusion 

Controlling systematic uncertainties in order to mea- 
sure sin 2 2#i3 with a sensitivity better than 0.01 at the 
90% confidence level is a challenge. A relative mea- 
surement using near and far antineutrino detectors can 
greatly reduce the detector-related systematic uncer- 
tainties. The Day a Bay experiment is designed with 
eight nearly identical antineutrino detectors. Multiple 
antineutrino detectors at each site enable side-by- side 
comparisons to estimate the relative uncertainties in de- 
tector efficiencies. The first two antineutrino detec- 
tors have been installed in EH1 and are operating with 
steady data- taking since September 23, 201 1. The anal- 
ysis of the first three months' data demonstrates a com- 
prehensive understanding of the detector responses and 
the relative detection efficiencies. The relative energy 
scale uncertainty is determined to be 0.4% and the rela- 
tive efficiency for the neutron energy selection is 0. 1 1 %. 
Combining with the precise target mass measurement 
and selection efficiencies, the relative neutrino detection 
efficiency uncertainty is 0.2% for the first two antineu- 
trino detectors, an improvement over the design value 
of 0.38% lit]. In this analysis, the expected ratio of the 
IBD rates in AD1 and AD2 is 0.981, compared to the 
measured ratio 0.987 ± 0.008(stat) ± 0.003(syst). 
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